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ABSTRACT 

Accreting neutron stars exhibit Type I X-ray bursts from both frequent hydro¬ 
gen/helium flashes as well as rare carbon flashes. The latter (superbursts) ignite in 
the ashes of the former. Hydrogen/helium bursts, however, are thought to produce 
insufficient carbon to power superbursts. Stable burning could create the required car¬ 
bon, but this was predicted to only occur at much larger accretion rates than where 
superbursts are observed. We present models of a new steady-state regime of stable 
hydrogen and helium burning that produces pure carbon ashes. Hot CNO burning 
of hydrogen heats the neutron star envelope and causes helium to burn before the 
conditions of a helium flash are reached. This takes place when the mass accretion 
rate is around 10% of the Eddington limit: close to the rate where most superbursts 
occur. We find that increased heating at the base of the envelope sustains steady-state 
burning by steepening the temperature profile, which increases the amount of helium 
that burns before a runaway can ensue. 

Key words: accretion, accretion disks — methods: numerical — nuclear reactions, 
nucleosynthesis, abundances — stars: neutron — X-rays: binaries — X-rays: bursts 


1 INTRODUCTION 

In low-mass X-ray binaries hosting a neutron star, hydrogen- 
and helium-rich material may be transferred from the com¬ 
panion star to the neutron star by Roche-lobe overflow. 
Compression due to a high surface gravity induces nuc¬ 
lear fusion, giving rise to a broad range of observed phe¬ 
nomena. Most prominent are Type I X-ray bursts with 
typical durations of 10-100s produced by runaway ther¬ 
monuclear burning (Grindlay et al. 1976; Woosley & Taam 
1976; Maraschi & Cavaliere 1977, see also Lewin et al. 1993; 
Strohmayer & Bildsten 2006). Stable burning is inferred 
from the absence of bursts at mass accretion rates above 
~ 10% - 30% of the Eddington limit (van Paradijs et al. 
1988; Cornelisse et al. 2003). In that case hydrogen and he¬ 
lium burn continuously in a steady state, without trigger¬ 
ing a thermonuclear runaway. All these burning regimes are 
reproduced by theory (Fujimoto et al. 1981), but there is a 
large mismatch in the conditions required for certain regimes 
to occur. For example, bursts are predicted to persist up to 
~ 100% of the Eddington limit (e.g., Heger et al. 2007a; 
Keek et al. 2009; Zamfir et al. 2014). 
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Many subsequent X-ray bursts or long periods of 
stable burning of hydrogen and helium burning leave be¬ 
hind an ashes layer that is rich in carbon. Runaway 
carbon fusion in these ashes is thought to power the 
rare “superbursts”: hours-long flashes from the same X- 
ray sources (Cornelisse et al. 2000; Cumming & Bildsten 
2001; Strohmayer & Brown 2002). All superbursting neut¬ 
ron stars also exhibit short bursts, but superburst obser¬ 
vations constrain the carbon mass fraction to be larger 
(~20%; Cumming et al. 2006) than what hydrogen/helium 
bursts are thought to produce (< 5%; e.g., Woosley et al. 
2004). Chemical separation could increase the carbon frac¬ 
tion (Medin & Cumming 2011), but would increase the re¬ 
currence times. These sources have relatively low burst 
rates (in’t Zand et al. 2004), suggesting that aside from 
bursts some stable burning occurs. This may produce larger 
amounts of carbon (e.g., Stevens et al. 2014), but no such 
stable burning mode was predicted by theory. 

In this Letter we investigate the ignition conditions of 
hydrogen and helium burning, and we identify a new stable 
burning regime at mass accretion rates bordering the burst 
regimes. Furthermore, time-dependent multi-zone models 
show that this regime is more prevalent when the base heat¬ 
ing of the neutron star envelope is strong. 
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2 NUMERICAL METHODS 

2.1 One-zone model 


We construct a simple one-zone numerical model to study 
the ignition conditions of hydrogen and helium burning as 
well as the stability of the burning processes in the neut¬ 
ron star envelope. The zone is given a temperature, T, and 
density, p, as well as mass fractions X, Y, Z, for *H, ^He, and 
the combined CNO isotopes, respectively. We use a solar ac¬ 
cretion composition of A = 0.73, Y = 0.25, and Z = 0.02, and 
we assume all accreted CNO is initially *^N. The Helmholtz 
equation of state and opacity routines (Timmes & Swesty 
2000; Timmes 2000) are employed to calculate the pressure, 
P, and opacity, k. The column depth is derived as y = P/g. We 
use a local gravitational acceleration g= 1.86 x lO'^cms^^, 
corresponding to a neutron star with a gravitational mass 
of IAMq and a radius of 10km in the Newtonian case or 
11.2km for General Relativity (e.g.. Keek & Heger 2011). 

Radiative cooling is implemented with a specific rate of 
(e.g., Fujimoto et al. 1981): 


acK 


( 1 ) 


We calculate the specific energy generation rate from nuc¬ 
lear burning, fnuc- For helium burning we use the Su rate 
from Fynbo et al. (2005). For hydrogen burning we em¬ 
ploy the hot /1-limited CNO cycle (/ICNO) for temperatures 
T > 0.8X 10^ K following Bildsten (1998), and at lower tem¬ 
peratures the CNO cycle using the *^N(p,7)*^0 rate as the 
slowest part of the cycle following Kippenhahn & Weigert 
(1994). Screening is implemented following Graboske et al. 
(1973). Finally, the base of the envelope is heated by 
pycnonuclear and electron-capture reactions in the neut¬ 
ron star crust (Haensel & Zdunik 1990, 2003; Gupta et al. 
2007; Schatz et al. 2014), or possibly by the dissipa¬ 
tion of rotational energy at the bottom of the ocean 
(Inogamov & Sunyaev 2010). The amount of heat generated 
per accreted nucleon, Q\,, is treated as a free parameter, and 
we include base heating as a specihc rate: q, = Qh^ily, where 
m is the specific mass accretion rate. We use a typical value 
for 2(, of O.lMeVu^* (e-g-j Gumming et al. 2006). The total 
heating rate is eheat = fnuc + eb- We assume the (small) con¬ 
tribution from compressional heating to be part of 2b (s-g-i 
Fujimoto et al. 1981). 

In the steady-state case, we determine m by equating 
the timescales for burning and accretion, Tnuc = Ucc, with 
Tacc = yhn-. 


m = ylTauQ- 


( 2 ) 


For helium burning Tnuc = FnucF/fnuc, where Enuc = 6.97 x 
lO^^ergg^* is the 3a energy yield per unit mass; for hydro¬ 
gen burning through the /ICNO cycle, Tnuc is set by the sum 
of the half lives of and *^0, t^ = 192.86s: Tnuc = Wr?’ 
where 14 accounts for the mass number of *^N, and 4 is the 
net number of protons that are burned in each cycle. We 
express m as a fraction of the Eddington limited value for 
solar composition: medd = 8.78 x 10^gs“*. 


dependent models of the neutron star envelope. The version 
of KEPLER and the setup of the models is similar to a pre¬ 
vious study by Keek et al. (2014) (see also Woosley et al. 
2004; Heger et al. 2007a, b). We refer to these studies for de¬ 
tails of the code, and here we summarize the main character¬ 
istics of our models. The code includes a large nuclear net¬ 
work using thermonuclear reaction rates from the REACLIB 
2.0 compilation (Gyburt et al. 2010). The network includes 
all CNO and jSCNO cycle reactions, the 3a process using 
the rate from Caughlan & Fowler (1988),^ as well as many 
other reactions such as break-out from the jSCNO cycle. We 
resolve the neutron star envelope in the radial direction. Ma¬ 
terial of solar composition is accreted onto a 2x 10^^g iron 
substrate, which acts as a thermal buffer into which heat 
from nuclear burning can dissipate. Crustal heating is im¬ 
plemented as an inflowing heat flux at the inner zone. The 
same accretion composition and neutron star parameters are 
used as for the one-zone model (Section 2.1). The results in 
this Letter are presented without GR corrections or redshifts 
(see Keek & Heger 2011). 


3 RESULTS 

3.1 Ignition conditions 

We use the one-zone model to reproduce the ignition 
conditions of hydrogen and helium burning derived by 
Fujimoto et al. (1981). Steady-state burning requires a bal¬ 
ance between heating and cooling: Cheat = fcool- For a range 
of temperatures, we determine the column depth where this 
condition is met for hydrogen and helium burning separately 
(“stable” in Fig. 1), and during stable burning the respective 
element is depleted at that location. A thermonuclear run¬ 
away ensues when the heating rate increases faster with T 
than the cooling rate. Both for hydrogen and helium burning 
we locate the boundary of the unstable region (“unstable” in 

Fig. 1): 

dcheat _ ^^cool 

dT ~ dT ^ ^ 

Material accreted onto a neutron star is compressed 
and heated over time, such that it moves to higher y and 
T in Fig. 1. The burning behaviour is determined by where 
the “stable” and “unstable” lines are reached (Fujimoto et al. 
1981; Fushiki & Lamb 1987; Bildsten 1998). Although CNO 
burning of hydrogen is unstable at low temperatures, over 
most of the considered temperature range hydrogen burn¬ 
ing is stable, because the /ICNO rate is independent of 
temperature. The stability of helium burning changes at 
T ^ 5x 10^K (we take the location to be where dTIdy = 0; 
e.g., Zamfir et al. 2014). Furthermore, as the stable hydro¬ 
gen line crosses the unstable helium line, helium flashes may 
occur either before hydrogen is depleted or after all hydro¬ 
gen has been burned to helium. Two helium ignition curves 
are, therefore, presented: one for the accreted helium frac¬ 
tion y = 0.25 and one where all hydrogen has burned to pro¬ 
duce Y = 0.98. The transition between pure-helium bursts 
and mixed hydrogen/helium bursts is at 7 == 2x 10^K. 


2.2 Multi-zone time-dependent models 

We employ the one-dimensional multi-zone hydrodynam¬ 
ics code KEPLER (Weaver et al. 1978) to create time- 


* Within the temperature range considered in our study, the 
Caughlan & Fowler (1988) and Fynbo et al. (2005) 3a rates differ 
at most 4%. 
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Figure 1. Column depth, y, as a function of temperature, T, for 
stable and unstable burning of H and He from a solar accretion 
composition, as well as a He mass fraction of 0.98 when all H 
has burned to He. Accreted material moves to larger y and T 
with time, and thus reaches any of the ignition curves. Starting 
at the stable H line, stable burning depletes He after Ay: the He 
depletion depth is the stable H line plus Ay. The thick solid line 
indicates where the new stable burning regime of H and He occurs 
when Ay <stable H line. The temperature profiles for two time- 
dependent models are shown: steady-state burning (upper line; 
2b = 0.75MeVu“*) and just prior to unstable He ignition (lower; 
2b = 0.1 MeV *). Their respective locations of H and He depletion 
are marked with ‘-h’. 

3.2 Helium depletion before runaway 

These ignition conditions are calculated assuming that hy¬ 
drogen and helium burning are independent. If only helium 
is accreted, the temperature increases with depth until igni¬ 
tion is reached. Because of the strong temperature depend¬ 
ence of the 3a rate, very little helium burns before the ig¬ 
nition depth. If, however, hydrogen burning takes place at 
shallow depths, the envelope is heated substantially, and he¬ 
lium burning is initiated already before its ignition curve 
is reached. Therefore, before helium burning dominates the 
energy generation rate, some helium already burns away. 

For each point on the stable hydrogen burning line 
(Fig. 1), we calculate the 3a burning rate assuming all hy¬ 
drogen has been converted to helium, and we determine the 
helium burning timescale, Tnuc (Section 2.1). Given the m 
associated with steady-state hydrogen burning (Eq. 2), we 
estimate the change in column during Tnuc as Ay = Tnuc'h. 
Here we assume that the temperature remains constant go¬ 
ing to larger values of y (e.g., Fujimoto et al. 1981), and we 
neglect the small increase in Cnuc from the increasing dens¬ 
ity. Where Ay is smaller than the depth of stable * H burning 
(Fig. 1), helium is depleted at a depth close to where hy¬ 
drogen burns: both hydrogen and helium burning proceed 
in a steady-state, and no thermonuclear runaway occurs. 
This is a new regime of stable nuclear burning. It happens 
only in a small temperature interval around T = 2.6 x 10* K, 
which is associated with a small range of accretion rates 



Figure 2. Helium mass fraction, Y, as a function of column depth, 
y, for two time-dependent models with base heating, 2b, in units 
of MeVu^'. They overlap at y < 6X 10^ gcm“^, but at larger depth 

Y drops faster for the hotter model. Dotted lines mark the condi¬ 
tions for unstable ignition for the two models. The colder model 
hits the unstable curve and ignites a burst, whereas the hotter 
model does not reach its unstable curve and has stable burning. 

0.08 < m/mEdd < 0.11 (thick solid line in Fig. 1). At lower T 
pure helium bursts ignite, whereas at higher T mixed hydro¬ 
gen/helium bursts occur. 

This is a simplistic approximation of time-dependent 
burning. Although it gives a good qualitative description 
of the behaviour, time-dependent models are required for 
accurate quantitative predictions as well as for a better un¬ 
derstanding of the new stable burning regime. 

3.3 Two time-dependent simulations 

For m = 0. IfitEdd and Qy, = O.lMeVu^*, where the one-zone 
model predicts the new stable regime, a KEPLER model ex¬ 
hibits mixed H/He bursts, indicating the new regime must 
occur at somewhat lower m for multi-zone models. There¬ 
fore, we create two KEPLER simulations with m = 0.02 mEdd^ 
one with the same 2b = 0. IMeVu^* as the one-zone models, 
and one with stronger base heating of 2b = 0.75MeV. As 
expected from Fig. 1, the first model exhibits pure helium 
bursts after hydrogen is depleted (similar to model Zm by 
Woosley et al. 2004). The second model, however, displays 
the new steady-state regime, even though m is much lower 
than where we expected it to occur. We compare the ignition 
conditions — which exhibit for 2b = 0.75MeVonly small 
shifts with respect to 2b = 0. IMeVu^* — to the T-profiles 
of the two models, where for the unstable model we select a 
time just prior to the runaway (Fig. 1). The locations where, 
respectively, the hydrogen and helium mass fractions are re¬ 
duced by an order of magnitude agree well with the one-zone 
estimates. 

Even though the two models have roughly similar tem¬ 
perature profiles, their burning behaviour is rather different. 
Up to the depth of hydrogen-depletion, the helium mass frac¬ 
tion is identical for both models (Fig. 2). For the steady-state 
model, the temperature profile is slightly steeper (Fig. 1), 
and the strong temperature dependence of the 3a reaction 
causes helium to burn faster than in the other model. As 

Y drops, the relevant unstable ignition curve moves to lar- 
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Table 1. Theoretical Nuclear Burning Regimes “ 
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Figure 3. Structure of the neutron star envelope as a function of 
column depth, y, from a time-dependent multi-zone simulation of 
steady-state burning. Top: specific energy generation rate, e„uci 
density, p, and temperature, T. Bottom: mass fractions for sev¬ 
eral isotopes. CNO is the sum for all C, N, and O isotopes in¬ 
cluding ’^C. At y <4x lO’gcm^^ CNO consists mainly of ’^O and 

* 50 . 


ger y. We create a series of ignition curves for a range of Y, 
and show their points of intersection with the T-profiles of 
the KEPLER models (Fig. 2). For the unstable model, Y is 
reduced by stable burning to 0.63 before it hits the instabil¬ 
ity line and ignites a burst (Fig. 2; see also Woosley et al. 
2004). For the steady-state model, however, stable helium 
burning proceeds fast enough to avoid its instability line, 
and a thermonuclear runaway is never reached. 

We study the steady-state model in more detail. 1.9 years 
of accretion and burning were simulated, during which time 
the burning layer was replenished 124 times. The helium 
mass fraction peaks at the depth where hydrogen is de¬ 
pleted, and rapidly decreases with increasing depth before 
the conditions for runaway 3a burning are reached (Fig. 3). 
The temperature in the helium burning region is below 
T = 2xlO*K, which means that the *^C(a,7)*®0 reaction 
is very slow compared to 3a, and a-capture is inefficient. 
Moreover, because almost all is produced at depths 
where hydrogen is depleted, *^C(p, 7 )*^N does not play a role 
either. The burning ashes, therefore, consist exclusively of 
98% in addition to the 2% accreted CNO mass fraction. 
The latter are predominantly and from the jSCNO 
cycle with a ratio of about 1 : 2. 


“ For solar accretion composition and Qb = O.lMeVu^'. 

*’ Peng et al. (2007), including sedimentation. 

Heger et al. (2007a). See also Keek et al. (2009); Zamfir et al. 
(2014); Keek et al. (2014). 

out (Fisker et al. 2007; Keek et al. 2014), rotational mixing 
(Piro & Bildsten 2007; Keek et al. 2009), and base heating 
(Keek et al. 2009; Zamfir et al. 2014). The small size of the 
m interval may explain why it has not been identified in 
observations. The mechanism that we describe to burn all 
helium stably occurs to a lesser extent in the regime of pure 
helium flashes, where part of the helium is burned before the 
thermonuclear runaway (Woosley et al. 2004). For accretion 
rates approaching the new stable regime, this part may be 
larger. 

In time-dependent multi-zone models we find that the 
same steady-state regime can take place at lower m, if the 
base heating is stronger. We used Qo = 0.75MeVu“*, which 
is in line with the amount of heat generated by electron cap¬ 
tures in the crust (Gupta et al. 2007). At lower m a higher 
value of 2b is actually expected, because of the temper¬ 
ature sensitivity of neutrino cooling (e.g., Gumming et al. 
2006). Furthermore, strong base heating may be needed 
to explain the observed recurrence times of superbursts 
(Keek & Heger 2011) and the ignition of superbursts in tran¬ 
sients (Keek et al. 2008; Altamirano et al. 2012). The tem¬ 
perature in the outer envelope is largely set by ySGNO burn¬ 
ing, whereas at the bottom of the envelope crustal heating 
sets an inflowing flux, such that a higher 2b leads to a higher 
temperature at the inner boundary. The temperature profile 
is, therefore, steeper for the model with a higher 2b' The 
boundary of the new steady-state regime is thus strongly 
dependent on 2b i and a large series of time-dependent or 
steady-state models are required to investigate it. 


4 DISCUSSION 

We have identified a new steady-state nuclear burning re¬ 
gime on accreting neutron stars. A range of burning be¬ 
haviour takes place as a function of and the new re¬ 
gime occurs in a small m interval close to 0.1/hEdd between 
the pure helium bursts and the mixed hydrogen/helium 
bursts. We summarize the different burning regimes in 
Table 1, although it should be understood that the bound¬ 
aries between regimes can be modified by, e.g., the metal- 
licity (Woosley et al. 2004), the details of GNO break- 


4.1 Superburst fuel 

Assuming all energy is provided by carbon burning, cool¬ 
ing models infer from superburst lightcurves a carbon mass 
fraction of 15%-26% (Gumming et al. 2006). Part of the 
energetics could originate from photodisintegration of heavy 
isotopes near the rp-process end-point (Schatz et al. 2003), 
but the composition of the ashes is likely dominated by iron- 
group nuclei (Woosley et al. 2004). Unstable burning pro¬ 
duces at most 5% (e.g., Woosley et al. 2004), whereas 

the new steady-state regime yields as much as 98%: neither 
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gives the amount inferred for superbursts. Superbursting 
sources have relatively high values of the ratio of the fluence 
in-between and during bursts (the so-called or-parameter; 
in’t Zand et al. 2004), which suggests that part of the ac¬ 
creted hydrogen and helium burns in a stable manner and 
part in bursts. To obtain 20% the steady-state mode 
must be active 16% of the time. As the steady-state occurs 
in a small m range, the 16% duty cycle may be provided by 
the intrinsic variations in m. Alternatively, there could be 
a hybrid between the bursting and stable modes, where a 
substantial part of ^He burns stably as in the steady-state 
regime, but an instability is triggered on a delayed times¬ 
cale. This regime would be different from the “delayed mixed 
bursts” proposed by Narayan & Heyl (2003). Further time- 
dependent simulations are needed to study the transition 
between the different regimes, and attention must be given 
to the possibility of some carbon destruction prior to super¬ 
burst ignition (e.g.. Keek et al. 2012). 


5 CONCLUSIONS 

We present models of a new steady-state nuclear-burning 
regime for hydrogen and helium mixtures accreted onto 
neutron stars. As hydrogen burns stably via the hot CNO 
cycle, the increased temperature induces Sor burning of he¬ 
lium. This mechanism fully depletes helium in steady-state 
at mass accretion rates of ~ 10% of the Eddington limit, 
bordered by burst regimes at lower and higher accretion 
rates. Enhanced base heating can steepen the temperat¬ 
ure profile and sustain the steady-state regime, extending 
the range of accretion rates where stable burning is active. 
Pure-carbon ashes are produced, which may power the rare 
energetic superbursts. In fact, more carbon is produced than 
needed to explain the fuel mixture of superbursts. Steady- 
state burning only needs to occur ~ 16% of the time. The 
new stable burning mode could contribute to the large values 
of the ff-parameter observed for superbursting sources. 

In future papers we will perform simulations to invest¬ 
igate the dependence of this burning mode on base heating, 
metallicity, and other parameters. Large burst catalogues 
such as MINBAR (Keek et al. 2010) will be employed to 
constrain the observational signature. 
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